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Fundamental properties of the Kosterlitz-Thouless-Berezinskii (KTB) transition which occurs 
in systems in the universality class of the two-dimensional X-Y model are reviewed here with an 
q \ emphasis on the real-space renormalization group analysis used to derive the recursion relations. A 

■ derivation of the unique temperature dependence of the correlation length of this system will be 

presented and the signatures of the KTB transition in electronic transport measurements will be 
, explained. Extensions of this model and technique to layered superconductors will also be presented. 

The size of the 3D critical region, the nature of the 3D crossover, and the effect of a finite electric 
current will be discussed. 
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The phase transition of the systems in the universality class of the two- dimensional (2D) X-Y model, known as the 
Kosterlitz-Thouless-Bcrczinskii (or some permutation of this) transition (Berezinskii 1971; Kosterlitz and Thoulcss 
1973; Kosterlitz 1974), is a fascinating one. And, even though the basic details of this transition were worked out in 
the early 1970's, interest in this transition remains keen with more than 200 articles published yearly on the topic. 
I | Indeed, one can find references to the "KTB" transition in the literature on topics ranging from layered materials 
(such as the high-temperature superconductors (HTSC's)) to supcrfluid films. Given its wide-spread importance, an 
understanding of the Kosterlitz-Thouless-Berezinskii transition and the derivation of its critical behavior properties is 
essential for a student of condensed matter and statistical physics. Renormalization group (RG) techniques (Wilson 
1971) are also an important tool for a theoretical physicist and can be used for a variety of problems in a number of 
ways. [For example, see the article by N. Andrei in this same issue.] It was originally applied to the 2D X-Y model 
by Kosterlitz (1974) to study its critical behavior. In this paper, we will review his real-space RG analysis of the 
2D X-Y model, describe the quintessential properties of the critical behavior, and discuss recent results on layered 
superconductors derived using the same mathematical formulation. 



II. THE 2D X-Y MODEL AND THE REAL SPACE RENORMALIZATION GROUP ANALYSIS 



Falling into the universality class of the 2D X-Y model are 2D superconductors and superfluids as well as 2D 
two-component magnetic systems. More generally, one can say that this universality class covers 2D systems with a 

S two-component order parameter. Systems in this universality class are believed to undergo the Kosterlitz-Thouless- 
Berezinskii phase transition. In this section, we will show how the 2D X-Y model can be mapped to the 2D Coulomb 
gas and then review Kosterlitz 's RG study of that system which illuminates many of the principal features of the 
critical behavior. (For pedagogical arguments of the phase transition see Kosterlitz et al., 1973 or Halperin (1979).) 
The RG study culminates in the 2D recursion relations which can be used to derive the temperature dependence of 
the correlation length. In the last subsection we will describe how KTB behavior can be detected in superconducting 
films via electronic transport measurements. 

One might guess that the physics of the 2D X-Y model is uninteresting since it has been shown that there is no 
long-range order at finite temperature in magnetic systems (Mermin and Wagner 1966) nor in superfluids (Hohenberg 
1967). Yet hints of novel behavior were provided by Stanley and Kaplan (Stanley and Kaplan 1966) who found a 
divergence of the susceptibility based on a high- temperature series expansion. Wegner (Wegner 1967) went on to 
show that the susceptibility was infinite at low temperatures. Further suggestions of a phase transition were provided 
by Berezinskii (Berezinskii 1971) who found a power- law decay of the spatial correlation function at low temperatures 
while others had ascertained the exponential decay at high temperatures. 

It was Kosterlitz and Thouless (Kosterlitz and Thouless 1973) who formulated the transition explicitly in terms of 
vortices and, to be more precise, as a vortex/anti- vortex pair unbinding transition. In this zero-field transition, vortices 
are created spontaneously in pairs at finite temperatures and have an energy which goes as the logarithm of their 
separation. As the temperature increases, the vortex pairs increase in size and in density. As this happens, screening 
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of their interactions becomes important and enables a phase transition at a temperature Tktb- At this temperature 
two things happen: vortex pairs start to unbind and spontaneous creation of single free vortices becomes possible. 
Being driven by topological excitations is one of the properties of this transition that makes it so unique. Another 
is its correlation length which diverges according to an exponential rather than a power-law at Tktb- As we will 
explain, it is this latter property that makes detection of this phase transition difficult to detect by thermodynamic 
measurements because there is no singularity in the free energy. Evidence for the transition can be accumulated 
however through electronic transport measurements for the case of superconductors. 



A. Mapping the 2D X-Y Model to the 2D Coulomb Gas 

The 2D X-Y model is defined by the Hamiltonian: 

H = - J s i ■ s j = - Js2 X! cos (^ ~ 03')) (!) 

<i,j> <i,j> 

where < i,j > are all nearest- neighbor spins and fa is the angle that the spin vector Sj makes with an arbitrary axis. 
In the low temperature and continuous limit, one has 

H = E + ^ J d 2 r\Vfa 2 , (2) 

where Eq is a constant. One can immediately see from Eq. (^|) the similarity to a Ginzburg-Landau free energy 
functional in which the amplitude variations are neglected. We now proceed to express Eq. (|J) in terms of the energy 
for a 2D Coulomb gas. This is accomplished (Halperin 1979) by first splitting the angle (j> into two parts, <f> = fa + <fi sw , 
corresponding to a rortex part and a spin wave part where 

V x Vfa w = 0, (3) 



V • Vfa = 0. (4) 

This substitution leads to three terms in Eq. (Q) including a cross term (J d 2 r\7(j) sw • V<^„) which can be shown to be 
equal to zero after an integration by parts leaving the vortex and spin wave contributions independent of one another: 

H = E a + ^- J d 2 r\Vfa\ 2 + J d 2 r\Vfa w \ 2 . (5) 

It is the vortex term that interests us because, as we will show, it contains the part that drives the phase transition. 

The next step in getting to the 2D Coulomb gas is to assume that vortices are present in the system which means 
that the line integral of the change in the field fa around any contour c is 2irN where N is any integer: 

■ dr = 2nN. (6) 

The simplest, non-trivial case is when the contour encloses exactly one vortex of vorticity one in which case N = 1. 
The right hand side (RHS) of Eq. (^) can be written in terms of a vortex density er(r) and the left hand side (LHS) 
can be rewritten using Stoke's theorem yielding 

z • (V x Vfa) = 27rcr(r). (7) 

By using vector identities and defining Vfa v = Vx Vfa, Eq. (0) becomes Poisson's equation: 

V 2 fa v = 2Tta(r). (8) 

Using a Green's function method to solve Eq. (^), one arrives at 

fa v (r 1 )= [ d 2 r 2 a(r 2 )l n | ri -r 2 |/£ (9) 



where £o is an infrared cutoff associated with the size of the core of the vortex. Finally, since Wfa ■ Vfa = Vfa v ■ "Vfa v 
and V • (4>' v \7 fa v ) = \7fa v ■ Vfa v + fa v V 2 fa v , one gets the energy for a configuration of vortices to be 
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d 2 r / d 2 r'a(r)a(r') In |r - r'j/Co 



d 2 m(r) / dV<7(r')ln(L/£ ), 



(10) 



where L is the size of the system. The first terms in Eq. ( |l0| ) contain the spin wave terms and the energy associated 
with the vortex core energies. The third term is an interaction term between the vortices and the last term which 
diverges as the size of the system ensures system neutrality. 

We can now neglect the spin-wave term, assume vortex neutrality, and write the partition function for the 2D 
Coulomb gas, 



N W V! J JDi JD 2 JD 2 n 



T2N 



x exp 



9 ^PiPjVQ r i ~ r j I) 



(11) 



where 2N is the total number of particles, N of which have a positive (negative) charge pi = +p (pi = —p), and 
are the coordinates of the ith charge. = fc^T where T is the temperature and fee is the Boltzman constant. 



y = cxp(/3/x)/r 2 is the fugacity, where \i = —E c and E c is the 
normal core in a superconducting fluid (in our case). V(R) = - 
notation ry = r,; 



'core energy," that is, the energy needed to create a 
ln(i?/r) is expressed in the units p 2 . Note that the 
• Tj will be used. The integrals are over an area Di which is all of the area except for disks di(j) of 



radius r around the charges j < i. The area Di is written, where A is the total area, 



Di = A- 



(12) 



B. RG Study of the Neutral 2D Coulomb Gas 

The renormalization group technique (Wilson 1971) is useful for studying second order phase transitions which are 
scale invariant at the transition temperature. RG incorporates the essential idea of scaling and extends it. Scaling says 
that near a second order phase transition, one length scale, the correlation length, dominates the critical behavior 
and all other length scales are thrown out. RG does not eliminate the short length scales altogether as scaling 
does, but incorporates their contribution into the renormalization of the system parameters. There are two steps 
in an RG analysis: the first is to course-grain the system which amounts to taking out the small-scale structure 
and incorporating it into a renormalization of the parameters of the system. This can be done both in real space 
and in momentum space. (Wilson and Kogut 1974). The second step is to rescale the system so that the infrared 
cutoff is restored to its original value. The connection to the physics of the system is achieved by studying how the 
parameters of the system change as one goes through the RG iterations. The relationship between the old values and 
the "renormalized" values are called recursion relations. For a thorough description of the renormalization group, 
the reader should consult any of the many textbooks on the subject (e.g. Binney, Dowrick, Fisher, and Newman 
1995; or Ma 1976). Below we sketch the details of the real space renormalization group study of the 2D Coulomb 
gas (Kosterlitz 1974). A manuscript containing the complete details of this RG analysis and the derivation of the 
recursion relations is available from the author. The reader is also encouraged to see any of the excellent reviews of 
the KTB transition (Halperin 1979; Minnhagen 1987; Suzuki 1979). 

The first step in the RG calculation is to integrate out small scale structure. In our case, this amounts to increasing 
the minimum separation between the vortices r to t + dr and determining the effect of the vortex pairs with a 
separation in that range on the interactions of other vortices. This is realized by first increasing slightly the allowed 
size of the disks around each vortex Di in Eq. ( [l2| ) from r to t + dr: 

f d 2 n = f dV< + E / «*V (1 3 ) 

When substituted into the integrals of Eq. (pi]), one gets, 
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(14) 



It is the last two integrals of the second term of the right hand side of Eq. (jTj) which concern us. By doing these 
integrals, the effect of the small pairs is integrated out. Physically, the integral over 8i(j) puts the charge i in an 
annulus around charge j to form a pair of smallest separation. The integral over D" then move this pair through all 
possible positions in the system. To do the integrals over and r 3 , all of the terms in the partition function that 
include i and j are isolated: 



D" 



d 2 rj / c^exp -@'^PjPkV(rjk) - py^PiPhVink) 



One can then assume that only vortices of opposite charge can form pairs so that pi = —pj, 

d^exp -j3pj /]Pk[- M r W T ) + ln(rf fc /r)] . 



d 2 r, 



D" 



Si(j) 



(15) 



(16) 



Integrating over r,, one gets = ry + r and, 



71 IT I d Tj 



d9 exp 



-#Pi£p*9 ln [(^ + 2r jk -T + T 2 )/r%] 



(17) 



One then makes the approximation that the density of vortices is small and consequently that rjk is large meaning 
that the probability that another vortex is close to the pair i-j is small. An expansion in r/rjk can then be performed, 



2nTdT J d i r j ( 1 + 



dO (f3pf 
2tt 2 



^2/PkPi 



rrjk cost/ rrjicosy 



r 



After the straightforward integration over 9, one has, 
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The integration over r,,- is also straightforward and one finds, 
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k=£l 



We can now write the second term of the RHS of Eq. (|lj) : 
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d n_i / d r i+ i. 



d fj_i / d r j+1 . 



2N-2 



d r 2 Arexp 



k^l:k.l^ij 



- ^ p fc p. In — 



2-KTdT 



2ttt 2 !3 2 P 2 x r fei 
A In — 



(18) 



(19) 



(20) 



(21) 



After a change of variables in Eq. (|2l]) (N' — N — 1), this equation can be added to the first term of the right hand 
side of Eq. <fiA 
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1 
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k^l 



T 



exp 



£/'-/'.. l:: 



(22) 
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Exponentiating the terms in the square bracket, 



Z = cxp[27ry 2 rdrA] ^ 



2N 



1 



U (Nl) 2 j D[ 



d 2 r x . 



D 



d 2 r 2N exp | - - ^2 PiP] 



('-M-fr)"?]) 



(23) 



Finally, we can rescale the lengths in Eq. ( |23| ) so that the limits in the integral are the same as those of the original 
partition function using r' = r/(l +dr/r), 



Z = exp[2wy 2 TdTA] ^ (y [l + dr/r(2 - /3p 2 / 2 ) 



27V I 



D 



N 

il 2 i]... I d 2 r 2 N exp | - - y^PiPj 



(AT!)2 



D 2 



l-(2^)^^).„^]} 



(24) 



One can see that this expression has the same form as the original except that the parameters have been renormalized: 

((3p 2 )' = Pp 2 (l - (2TTyT 2 ) 2 Pp 2 dr/2r), (25) 



y' = y(l + (2 - (3p 2 /2)dT/ T ) 
These equations can be put in differential form, 



d{Pp 2 ) 
(PP 2 ) 2 



-d{f3p 



2\-l 



(2TTyT 2 ) 2 dr 
27 



a y 131)2 / 4 i.w / 
dy = — { W 2 ~ )lT - 

and the definitions e = ln(r/£o), y' — 2nyr 2 , and K — (3p 2 can be made yielding, 

dK/de = -y 2 K 2 /2, 

dy/de = y(4-K)/2, 

where we have dropped the primes. These equations are frequently written in terms of x = A/ K — 1: 

dx/de — 2y 2 , 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



dy/de = 2yx, (32) 
where we have linearized the equations in x. These equations can be integrated to find x(e) (or K(e)) and y(e). 



C. Analysis of the Recursion Relations 

The recursion relations Eqs. (|3l])-(|3^) have a fixed line at y — along which the parameters of the system do not 
change. But it is the fixed point [x = 0,y = 0] that is of more interest to us as we shall see. How these parameters 
change as one goes through the RG iterations are represented as "flows" in phase space. By inspecting Eq. (|3l]), one 
sees that x can only increase under a RG iteration but that the increments can become infinitesimally small as y goes 
to zero. Eq. (|3|) on the other hand tells us that y increases when x is positive and decreases when x is negative. 
This behavior is reflected in Fig. |l|a where the RG flows are plotted in x-y space for various initial values of x and 
y. One can see that the lines end on the y axis, tend toward [x = oo,y = oo], or, in one case, flow to [x = 0,y = 0]. 
This reflects the tendency of the RG iterations to drive the system to its "pure" limits: the low-temperature limit 
of no vortices (y = 0) or the high- temperature limit of a high vortex density. The last case where the flows go to 
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the fixed point represents a sepatrix or, in RG parlance, a critical surface. This surface is associated with the critical 
temperature (where the system is scale- invariant) since it does not flow to either the low or high temperature limits. 
Plotted in Fig. [l]b are the RG flows for Eqs. (|2^)-(|3^) with initial values of K and y corresponding to those used in 
Fig. 0a. The same behavior is reflected. 

y 2 - x 2 = c, (33) 



The first integral of Eqs. M-(B2ft is 



where c is a constant. Eq. (33) is the equation for a hyperbola and the flows of Fig. [l] represent a family of hyperbolae 
corresponding to different values of c. Based on our discussion in the above paragraph, it is natural to make the 
assumption, 

c cx (T - T KTB ) / T KTB ■ (34) 

For c > 0, the system is at a temperature above the critical one and the flows go to the high-temperature limit. For 
c = 0, T = Tktb and the RG flow goes to the origin. For c < 0, one is in the low-temperature limit. 

With Eq. (p3|) and Eq. (|34| ), one can derive the temperature dependence for the vortex correlation length. Substi- 
tuting Eq. (|33|) into Eq. Q31|), one has dx/de — 2(x 2 + c). Integrating, one gets 



1 x 
— arctan — — 



2e, (35) 



where X\ is the bare or initial value of x and x rnax is the final value. To be more precise, one integrates the recursion 
relations until there are no longer valid which is when the vortex density becomes large (i.e., when y becomes of 0(1)). 
The value of e at that point is e max . Because the correlation length £ is dominating the critical behavior and because 
the only length scale in Eq. ( p5| ) is r = £,o exp[e max ] , one must make the association £ = r. Then, for small c, one 
finds 



£(T) cx exp[y/b/(T-T K TB)]. (36) 

This result is valid for T > Tktb and for T < Tktb one can show that the correlation length is infinite (Kosterlitz 
1974), which agrees with the result mentioned earlier that the susceptibility is infinite at small temperatures. This 
unique temperature dependence of the correlation length is a distinguishing property of the KTB transition. Typically, 
one has a power law divergence of the correlation length. The correlation length can be thought of as a screening 
length for the Coulombic interaction. For separations less than it, the interaction is logarithmic but approaches a 
constant at larger separations due to the screening of free vortices. It can also be interpreted as the average distance 
between free vortices. 

Further insights into the KTB transition can be gained by looking at the length scale dependence of K and y as 
derived from Eqs. (p9|)-(pc|). In Fig. ^, we have plotted K(e) and y(e) for the same initial values of K and y as those 
used in Fig. [l]b. From the length scale dependence of y(e), one can see that there are two types of behavior as in 
Fig. [l]b, one corresponding to T > Tktb and the other to T < Tktb- In the latter, y(e) falls monotonically to zero 
meaning that the density of vortices with a separation R — £olne becomes negligible at large separations. Above 
Tktb-, this quantity initially decreases meaning again that the density of vortices falls with increasing separation, but 
then starts to increase. This increase marks the onset of unbinding of vortex pairs at the larger separations and the 
spontaneous creation of free vortices. This behavior is reflected in K(e) depicted in Fig. 2b. For temperatures less 
than Tktb, K decreases due to the screening of vortex pairs but saturates at a finite value. The length scale at which 
it saturates corresponds to the length beyond which vortex pairs with the corresponding separation are negligible. 
The effective interaction is proportional to K(R) ln(i?/£o which remains largely logarithmic. This is in contrast to the 
T > Tktb case where one can see that K does not saturate but decreases to zero. This is because vortices are present 
at all length scales. In fact, it is not just K that changes but also the logarithmic dependence of the interaction energy 
due to free vortices. For further information, see the discussion by Minnhagen (1987). 

D. Observing the KTB transition in superconducting films 

An intriguing consequence of the KTB temperature dependence is that there is no singularity in the free energy of 
this system. Using a scaling analysis, one finds that the singular part of the free energy density scales as 

/ cc r 2 (37) 



G 



which is analytic as T — ► Tktb- The specific heat and the magnetization, both thermodynamic derivatives of the 
free energy, are also analytic at the critical temperature leaving detection of this transition difficult. 

In superconducting films, it is primarily through electronic transport measurements that KTB critical behavior is 
verified. Here we will briefly sketch the arguments leading to the current-voltage characteristics and the temperature 
dependence of the resistance R(T). In the presence of a current, a force (whose magnitude is linear with current 
and whose sign is opposite for vortices and anti- vortices,) is exerted on a vortex in a direction perpendicular to the 
current. Free vortices can therefore dissipate energy. The resistance in the limit of zero current should be proportional 
to the density of free vortices which in turn is inversely proportional to the square of the correlation length: 

R(T) (x n F (x C 2 ■ (38) 

Minnhagen (1987) has generalized this formula to include the effects of the underlying superfluid to obtain R(T) oc 
exp[-\/6(T c o — T)/(T — Tktb)] where T c q is the mean-field transition temperature. 

The effect of a current on a pair is to pull them apart with a force which is independent of their separation. Thus, 
the bare energy of a vortex pair in the presence of a current is 

E(R) = 2E C + p 2 logOR/G)) ~ IR/to (39) 

where we have expressed the current in units of energy. At large R, the current repulsion will always dominate and 
E(R) will have a peak at R c = £,op 2 /I. Thermal energy will activate hops over this barrier providing a mechanism 
for vortex pair unbinding with a rate (Huberman, Myerson, and Doniach 1978) 

r oc exp[-E B /k B T], (40) 

where Eb = 2E C — p 2 +p 2 \og(p 2 / 1) is the height of the barrier. The density of free vortices in this case can be found 
using the kinetic equation, 

dn F /dt = T -n%, (41) 

where the second term on the right hand side account for vortex/anti- vortex pairing or annihilation. In the steady 
state, the RHS is zero so that Up = T 1 / 2 . This results in a non- linear I-V relation: 

VocI a ^+\ (42) 

where, neglecting renormalization effects, a(T) — p 2 /2k B T. In the limit of weak currents, a(T) is decreasing linearly 
to a value, cx{Tktb) = 2, as one increases the temperature to Tktb at which point it jumps discontinuously to zero. 

The I-V signatures and the temperature dependence of R{T) that characterize KTB behavior, Eqs. ( |3^ ) and 
have been observed in superconducting films (Minnhagen 1987) providing strong evidence for the influence of vortex 
pair unbinding in these systems. They have also been observed in layered superconductors with a generalized form 
of Eq. ( f42| ) (Jensen and Minnhagen 1991) indicating that vortices also play a major role there. We now proceed to 
discuss the layered system. 



III. CRITICAL BEHAVIOR OF VORTICES IN LAYERED SYSTEMS 

In this Section, the model and techniques of Section || are generalized to a weakly coupled layered system. The 
various properties of the critical behavior of vortices in layered systems deduced using the methods reviewed above 
will be emphasized but it should be noted that many authors (Chattopadhyay and Shenoy 1994; Fischer 1993; Friesen 
1995, Horovitz 1991; Korshunov; 1990) have made significant contributions to the field using other approaches. We 
will begin with the expected behavior of a layered system and the basic model used for such systems and discuss 
how the bare interactions of the vortices are modified in a layered system. The recursion relations for the interacting 
layered vortex gas derived from a RG analysis will then be given and an analysis of these recursion relations will be 
used to study the three-dimensional to two-dimensional crossover. Finally, the effect of a uniform electric current I 
on the system and the I-T phase diagram will be considered. 

An intuitive understanding of the conventional wisdom (Leggett 1989) of the critical behavior of vortices in layered 
systems can be achieved by considering the correlation length which diverges at the critical temperature T c . When 
the correlation length is smaller than the distance d separating the layers, 2D behavior of the vortices is expected. 
When the correlation length becomes larger than the interlayer separation, the behavior should be 3D. Since the 
latter condition is met as the correlation length diverges near the transition, 3D behavior is expected in a small 
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temperature window around T c . Outside of this window, the layers are expected to be uncoupled. One can use 
£ c = A exp[y/b/(T — Tktb)] as an estimate for the correlation length in the c-direction (direction normal to the 
layers) where A is the ratio of the Josephson coupling energy to the in-plane coupling energy. Using the criterion 
£ c = d for the dimensional crossover, one obtains for the size of the 3D critical region t^jj a functional dependence, 
t~3d oc (lnA)~ 2 . We will see below that this argument breaks down above T c where vortex fluctuations proliferate. 

The most common model for a stack of weakly coupled superconducting layers is the Lawrence-Doniach model 
(Lawrence and Doniach, 1971) in which the coupling between the layers is taken to be the same as that between two 
superconductors separated by an insulator: Josephson coupling. The effect of the Josephson coupling on the vortex 
interactions is to strengthen the intralayer vortex interaction and to introduce interlayer vortex interactions. These 
interactions have been derived with the intralayer interaction V(R, 0) (Cataudella and Minnhagen 1990) being 



V(R,0) = 



- ln(i?/r) + (XR 2 /At 2 ) ln(Ai? 2 /r 2 ), r < R < R x 
~(nVXR/Ty/2), i?>i? A 



f -Ai? 2 /r 2 ln(Ai? 2 /r 2 ), R < R x 
\ VXR/t, R > R x 



and the interlayer interaction V(R, 1) (Bulaevski, Meshkov, and Feinberg, 1991) being 

V(R,1) oc 

where Rx = £,o/V\- For generality, one can introduce the effect of the current through the term, 

H = J2p i T i -3 ) (43) 

i 

where is the in-plane Cartesian coordinate of the ith charge, Pi = ±y/ 'd/2$o/27rA is its charge, and J is related to 
the current / by J = ^^dl /pea where d is the distance between the layers, <£>o is the superconducting flux quantum, a is 
the characteristic cross-sectional area that relates the current to the current density I /a, A is the London penetration 
depth, and c is the speed of light. 

The partition function for the layered vortex gas (Pierson 1994; Pierson 1995b) in the presence of a current is, 



Z 



e^vWe/ ^e/ ^ 2 -e/ ^ 

N h Jd i l 2 Jd * 1 2n Jd ^ 

x exp [ - f J2 Pip i V U r ' J ~ v i l> li ~ l i) + Pj2 PiTi ' J ] ' ( 44 ) 

where is the index of the layer in which the ith particle is located. Interactions between vortices separated by more 
than one layer are neglected. A renormalization group analysis has been carried out on this partition function and 
the following recursion relations were found, (Pierson 1995a; Pierson 1995b) 

dx/cfe = 2y 2 [l - A/16 + J 2 /(l + x)], (45) 



dy/de = 2[x+ (1/2) A In A]y/(1 + x) + Jy/VTT^, (46) 
d\/de = 2A{1 - 4y 2 [l + J 2 /(l + x)]/(l + x)}, (47) 

dJ/de = J, (48) 

where a factor of \fj3r is absorbed into J. Note that these recursion relations reduce to Eqs. ([3l])-([32]) for J = and 
A = 0. One can see that in Eqs. (^)-(^) the current tends to counteract the effect of the Josephson coupling. This is 
expected because the current introduces a repulsion between vortices of opposite sign while the Josephson interaction 
strengthens the interaction. 

Considering first the zero current recursion relations, the critical behavior of vortices in a layered system ca n be 



addressed. This is done by examining the correlation length which is defined in the same way as in Section II C : 
£(T) — £o exp[e m ax}- The correlation length will diverge like the 2D dependence sufficiently far above the transition 
temperature T c but will cross over to a power-law temperature dependence closer to the transition temperature. T c 
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is easily derived from this quantity since it peaks at this temperature. By varying A one can determine the functional 
dependence of T c on the interlayer coupling to find (Pierson 1995a) 

T c <x l/(ln(A)) 2 (49) 

in agreement with the results of Hikami and Tsuneto (1980). By comparing the correlation length to the 2D correlation 
length, insights into the size of the 3D critical region can be found. One will find that sufficiently far from T c the two 
correlation lengths are nearly identical which is interpreted to mean that the layered system behaves two-dimensionally 
in that region. Closer to T c however there is a slow bifurcation which is taken to indicate crossover to the 3D region 
and the size of this region can be studied as a function of A. ft is convenient to divide the 3D temperature window 
into two parts, one above T c : Tq D , and one below T c : t^ d . It is found (Pierson 1995a) that while t^ d (A) has the same 
logarithmic dependence on A as previously found for t$d, t£ d differs markedly: t^ d oc A 1 / 4 . A power-law dependence 
on A was also found subsequently by Friesen (1995): t£ d oc A 3 / 4 . The anisotropy of the 3D behavior around T c is 
dramatic as evidenced not only by the different dependences on A but also by the magnitude difference above and 
below the transition temperature (Pierson 1995a). This anisotropic behavior is attributed to the importance of the 
vortex fluctuations above T c which screen out both the interlayer and intralayer coupling. 

How the system crosses over the 3D to 2D behavior can also be addressed (Pierson 1995c) by studying Eq. d47|), 
the recursion relation for A. Below T c where y is small, A grows to the isotropic value. This means that as the size 
of the vortex pairs becomes larger, 3D effects become more important. Above T c , A grows initially but then starts 
to decrease at larger length scales due to vortex screening represented in the term proportional to y 2 . The length 
scale at which A(e) peaks is the interlayer screening length I^d/id an d it is found to get larger as one lowers the 
temperature towards T c presumably diverging there as the intralayer screening length (i.e. the correlation length) 
does. The implication of this for the nature of 3D to 2D crossover is that the layer decoupling occurs first at the 
largest length scales just above T c and then proceeds to smaller length scales as the temperature is increased. As a 
consequence, the interaction between two vortices in neighboring layers separated by a distance larger than 1 3 d/2D 
will be strongly screened giving the effect of decoupled layers. The linear (or 3D-like) interaction between two vortices 
in the same layers separated by a distance larger than interlayer screening length will also be screened resulting in the 
same effect. The screening of 3D effects at large length scales has interesting effects on the effective dimensionality of 
vortices where the 3D effects are felt primarily at those same lengths (Pierson 1995c). 

An analysis of the finite current recursion relations has also shown to be fruitful (Pierson 1995b). By comparing 
the correlation length to the zero current correlation length and the 2D correlation length, one finds that there are 
three characteristic currents. The smallest ll(T) is the current that marks the onset of resistance. It corresponds 
to the minimum current needed to overcome the Josephson attraction to allow vortex pair unbinding. The next 
characteristic current I C {T) is where the phase transition occurs and above which free vortices can be spontaneously 
created. The largest characteristic current I%(T) is the 3D/2D crossover current. The temperature dependencies of 
these characteristic currents have been examined and are found to be predominantly linear. While I^iT) is linear at 
larger currents, it crosses over to join I C (T) at small currents converging to the same point at T c . These three currents 
map out a phase diagram in I-T space. The properties of electronic transport properties in each region of the I-T 
phase diagram are discussed by the author (Pierson 1995b). It should be stressed that I^(T) and I 2 {T) are crossover 
currents while I C (T) is a phase transition. 

IV. SUMMARY 

In this paper, we have reviewed the critical behavior of the 2D X-Y model universality class with an emphasis 
on the derivation of the recursion relations through a renormalization group analysis. We reviewed the derivation of 
the temperature dependence of the correlation length and the methods of verifying the transition in superconductors. 
Finally those methods were generalized to a layered superconductor system where the vortex interactions are modified 
due to the Josephson coupling between the layers. The power and usefulness of the real space RG method become 
readily apparent here as evidenced by the calculation of a number of important properties of the critical behavior 
including the dependence of the transition temperature and the size of the 3D critical region on the strength of the 
interlayer coupling. 

ACKNOWLEDGMENTS 

This work was supported by the Office of Naval Research. 



9 



* Electronic address: pierson@wpi.edu. 
V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). 

J. J. Binney, N. J. Dowrick, A. J. Fisher, and M. E. J. Newman, The Theory of Critical Phenomena: An Introduction to the 

Renormalization Group. (Oxford University Press, New York, 1995). 

L. N. Bulaevski, S. V. Meshkov, and D. Feinberg, Phys. Rev. B 43, 3728 (1991). 

V. Cataudella and P. Minnhagen, Physica C 166, 442 (1990). 

B. Chattopadhyay and S. R. Shenoy, Phys. Rev. Lett. 72, 400, (1994). 

K. H. Fischer, Physica (Amsterdam) 210C, 179 (1993). 

M. Friesen, Phys. Rev. B 51, 12786 (1995). 

B. I. Halperin, in Physics of Low- Dimensional Systems, Proceedings of the Kyoto Summer Institute, 1979, edited by Y. Na- 
gaoka and S. Hikami (Publication Office, Prog. Theor. Phys., Kyoto, 1979) p. 53. 
S. Hikami and T. Tsuneto, Prog. Theor. Phys. 63, 387 (1980). 
P. C. Hohenberg Physical Review 158, 383 (1967). 
B. Horovitz, Phys. Rev. Lett. 67, 378 (1991). 

B. A. Huberman, R. J. Myerson, and S. Doniach, Phys. Rev. Lett. 40, 780 (1978). 

H. J. Jensen and P. Minnhagen, Phys. Rev. Lett. 66, 1630 (1991). 

S. E. Korshunov, Europhys. Lett. 11, 757 (1990). 

J. M. Kosterlitz, J. Phys. C 7, 1046 (1974). 

J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 

W. E. Lawrence and S. Doniach, in Proceedings of the International Conference on Low-Temperature Physics, Kyoto, 1970, 
edited by E. Kanda (Keigaku, Tokyo, 1971), p. 361. 

A. J. Leggett, in The New Physics, edited by P. Davies (Cambridge University Press, Cambridge, England, 1989). 

S.-K. Ma, Modern Theory of Critical Phenomena, (Benjamin/Cummings, Reading, MA 1976). 

N. D. Mermin and H. Wagner [PRL 17, 1133 (1966). 

P. Minnhagen, Rev. Mod. Phys 59, 1001 (1987). 

S. W. Pierson, Phys. Rev. Lett. 73, 2496 (1994). 

S. W. Pierson, Phys. Rev. B 51, 6663 (1995). 

S. W. Pierson, Phys. Rev. Lett. 74, 2359 (1995). 

S. W. Pierson, Phys. Rev. Lett. 75, 4674 (1995). 

H. E. Stanley and T. A. Kaplan PRL 17, 913 (1966). 

M. Susuki, in Physics of Low-Dimensional Systems, Proceedings of the Kyoto Summer Institute, 1979, edited by Y. Nagaoka 

and S. Hikami (Publication Office, Prog. Theor. Phys., Kyoto, 1979) p. 37. 

F. Wegner, Zeitschrift fur Physik 206, 465 (1967). 

K. G. Wilson, Phys. Rev. B 4, 3174 (1971). 

K. G. Wilson and J. Kogut, Phys. Rep. C 12, 75 (1974). 

FIG. 1. (a) The RG flows for the recursion relations Eqs. (pil] )-(|32"|) for various initial values of x and y [xi — —0.5, 
t/i = 0.44, 0.46j__- 0.56]; and (b) The RG flows for Eqs. @-(|(|for the initial values: K t = 8.0, yi = 0.56, 0.58, •■■ 0.68. 
Because Eqs. (|l|)-(|3|) have been linearized around x = 0, the initial values in (a) are not the same as those in (b). In both 
figures, there are two classes of flows, one for which y goes to zero and one where y becomes large, separated by a curve 
corresponding to T = Tktb- 

FIG. 2. The length scale dependence of (a) y and (b) K plotted versus e = ln(r/^o)- In the low temperature phase, K (e) 
saturates to a finite value at a length corresponding to the separation beyond which the vortex pair density becomes negligible. 
Above Tktb where the number of vortices increases at larger lengths, K(e) is renormalized towards zero. 
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